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Abstract. We discuss the recent progress in calculating the properties of 'hybrid stars' (stellar 
objects similar to neutron stars, classified by the incorporation of non-nucleonic degrees of freedom, 
including but not limited to hyperons and/or a quark-matter core) using the octet-baryon Quark- 
Meson Coupling (QMC) model. The version of QMC used is a recent improvement which includes 
the in-medium modification of the quark-quark hyperfine interaction. 
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1. INTRODUCTION 

The study of the QCD phase diagram is of great interest to the scientific community. 
Lattice QCD produces simulations which provide an insight into the {T > 0,jU ~ 0} 
region of the phase diagram, whilst properties of QCD at finite chemical potential can 
be calculated using various models for dense hadronic matter (at T > 0, /i ^ 0). In order 
to qualify the success of these models we require a method to test the predictions via 
experiment and observation. For finite chemical potentials, this requires tests at various 
relevant density scales. For low-density predictions we can compare calculations to 
experiments involving finite nuclei, but no experiment (to date) can probe the entire 
range of extreme densities believed to exist within neutron stars ^ We therefore turn our 
attention to proxy measurements which may support or refute our predictions, such as 
the masses and radii of pulsars. 

We focus our attention on a particular model — the Quark-Meson Coupling model 
(QMC), which we shall describe in Section 2 — that has seen much success at each of 
these density scales, and focus in particular on the high-density Equation of State (EoS) 
for infinite matter (see Section 3). In order to examine the effects of quark degrees of 
freedom at extremely high densities we model a phase transition between this hadronic 
EoS and a quark-matter EoS as described by several models (see Section 4). We examine 
the predictions for neutron star properties that arise from calculations involving the QMC 
EoS and the effect that the phase transition has on these properties in Section 5. 



In Ref. [1] the authors note that heavy-ion collisions may be able to provide insight into the properties 
of matter up the density range of 2-3po. 



2. QMC 



The Quark-Meson Coupling (QMC) model [2] describes baryons as 'bags' of three 
confined quarks for which the energy-density is greater than that of the surrounding 
non-perturbative vacuum by 5 ~ (180 MeV)4. These quarks are immersed in a mean- 
field of mesons which contribute to attractive scalar and repulsive vector potentials, 
which describe interactions in a similar way to Quantum Hadrodynamics (QHD). This 
model has seen much success in predicting the properties of finite nuclei, in particular 
hypernuclei (containing strange quarks) with the recent addition of a A-L hyperfine 
splitting contribution [2]. 

We model baryons influenced by scalar-isoscalar o, vector-isoscalar (O, and vector- 
isovector p mesons, noting that the mean-field pseudoscalar-isovector n field contribu- 
tion is zero due to parity considerations. We additionally include leptons (£ e jU~}) 
interacting with baryons in beta-equilibrium in our calculations. 

In QMC, the effective mass for the baryons has a quadratic form — incorporating the 
scalar polarizability d which self-consistently allows feedback of the scalar field — given 



where the factors Wg and Wg (the origins and values of which are discussed in further 
detail in Ref. [2]) account for the SU(6)-style coupling ratios for the various octet 
baryons B G {A^^^^*^, AjZ^'^jS^'*^}. The inclusion of hyperons into these calculations 
does not result in negative effective baryon masses — the occurrence of which is a 
shortcoming of QHD attributed to the linear effective mass, and the reason that we do 
not model hyperons in QHD — since the scalar potential responds to the strength of the 
scalar field in a way that prevents this. Indeed, the effective masses of all the baryons 
remain positive for any value of (a) in QMC. 

The couplings of the nucleons to the mesons in QMC (gN(y:8Nco-igp) are fitted such 
that the saturation properties of nuclear matter are reproduced in isospin symmetric 
nuclear matter. The nucleon-meson couplings and quantities that we reproduce are given 
in Table 1. The couplings of the hyperons to the isoscalar mesons are related to these 
nucleon-meson couplings via SU(6) spin-flavor relations, and while the couplings of the 
hyperons to the isovector meson are unified into a single value, each vertex term includes 
an isospin factor which splits the vector potentials within an isospin group. 



by 




(1) 



TABLE 1. Meson-nucleon couplings and reproduced saturation proper- 
ties of nuclear matter: energy per baryon; baryon density; and symmetry 
energy at saturation. The couplings are fitted to best reproduce the experi- 
mental values. 



8aN 8 (ON gp I / Po-Y.bPbMb) Po <H 



8.278 8.417 8.333 -15.86 MeV 0.16 fm"^ 32.5 MeV 



The hyperon Fermi momenta (and hence, densities) are determined via relations 
between the baryon chemical potentials, according to 



Hi = BiHn - Qifie = ^kl + {Mi + ^)^ + i:f, (2) 

where 5, = +1 \/ i E B and Qi are respectively the baryon and electric charges (the 
latter in units of the proton charge) associated with the neutron and electron independent 

s 

chemical potentials jU„ and ju^; and L^' are the scalar and temporal components of the 
baryon self -energy. 

In QHD, the baryon self-energies at Hartree level are calculated by considering the 
tadpole diagrams as modifications to the baryon propagator. Following the derivation of 
Ref. [3], the scalar component of the baryon self -energy is given by 

^BQis^ = -SaB{o) = -8aBl,-^—^^^ J ®{kF,,-\k\)^^d k, (3) 

1 /2 

where E^{k) = {k^ + {M%Y) ' , and where 7/ = i V i e 5 is the spin of baryon /. In 
QMC, the scalar self -energy is modified to obtain Eq. (1), which requires care to be 
taken with the additional self-consistency. 

Using these relations and fitting the couplings to saturation properties, the only inde- 
pendent quantities that still require definition are the proton and neutron Fermi momenta. 
These are constrained by requiring that the conserved total baryon density is equal to a 
selected value, and that the total electric charge of the system is zero. These conditions 
are given by 

Ptotai = Ep5, = £ QbPb + £ QePe- (4) 

B B £ 

3. INFINITE BARYONIC MATTER 

Once the required self-consistent equations of the previous section have been satisfied, 
we can calculate the properties of dense matter. The species fractions Yi of various 
particles at a wide range of densities can be calculated once all of the relevant Fermi 
momenta are determined, given that 

Yi = pi/p, pi=^^^0lj@{kF-\k\)d\=^^. (5) 

The energy-density and pressure in QMC are analytically identical to those of QHD (as 
given in Ref. [3]), noting of course that now the effective mass has the quadratic form of 
Eq. (1). 

We have assumed up to this point that at extremely high densities (p > 3po) hyperons 
will remain the correct degrees of freedom, however it is possible that at some density 
the concept of distinct baryons becomes inaccurate, and that the component quarks can 
become deconfined. In this case, we must consider the possibility of a phase transition 
to quark matter. 



4. PHASE TRANSITIONS 



The Gibbs conditions for a first-order phase transition between two phases (/, j) are that 
at the phase transition point, the following equiUbria are achieved: 

Thermal Equilibrium: The total temperature of each phase must be equal, thus 7] = Tj, 
Mechanical Equilibrium: The total pressure of each phase must be equal, thus Pi = Pj, 
Chemical Equilibrium: The independent chemical potentials of each phase must be 

equal, thus {iia}i = {l^ajj- 

It is possible that there exists a range of densities at which these conditions are met (and 
thus a mixed-phase exists), outside of which the phase with the greater pressure will be 
energetically favoured. Further details can be found in Refs. [4, 5] 

We calculate the baryon phase EoS at T = at increasing densities (using the inde- 
pendent chemical potentials as inputs to the quark-matter EoS) until we reach a point 
at which the pressure in each phase is equal (if it exists). This signals the onset of a 
mixed phase which we parameterize with the quark fraction < < 1 up to the phase 
transition to pure quark- matter at = 1. We then continue to calculate the quark-matter 
EoS with increasing density. 

We use two quark-matter models; the MIT bag model [6] and the Nambu-Jona- 
Lasinio (NJL) model [7]. The former models constant current-quark masses, while the 
latter models the effect of dynamic chiral symmetry breaking, producing constituent- 
quark masses at low densities, and current-quark masses at large densities. 

The species fractions for a hybrid EoS involving a hyperonic QMC baryon phase and 
an MIT bag model (using 

51/4 

= 180 MeV) quark phase are shown in Fig. 1 where we 
note the presence of baryon-, mixed-, and quark-phases. 




FIGURE 1. Species fractions Yi for hyperonic QMC with a phase transition to a MIT bag model quark 
matter via a mixed phase. We note that with these parameters, the transition to a mixed phase occurs below 
saturation density, po = 0.16 fm^^. 



If we calculate the hybrid EoS using hyperonic QMC and NIL we find that the quark 
pressure remains lower than the baryon pressure for any density. This is a result of 
the large quark masses at low density and the relation between the pressure and Fermi 
momentum (via /i^). As a consequence, we find that no phase transition from hyperonic 
QMC to NIL quark-matter is possible for any reasonable value of B. We do however 
find that the transition is possible from nucleonic QMC to NJL quark-matter, in which 
case the baryon phase is not significantly softened by the inclusion of hyperons. 



In order to test these models against observations, we calculate properties of hybrid stars 
such as the total stellar mass and radius. The masses of these objects are comparatively 
easy to measure in the form of (particularly, binary) pulsars. The radius however has 
proven to be a greater challenge to measure. 

In order to calculate the predicted properties of hybrid stars we use the EoS (which 
may include a mixed- or quark-phase) as input to the Tolman-Oppenheimer-Volkoff 
(TOY) equation 



where S' and P are the energy density and pressure in the EoS, respectively. 

The predicted values of mass and radius for QMC both with and without a phase 
transition to MIT bag model quark-matter are shown in Fig. 2. We note that the inclusion 
of a phase transition lowers the stellar mass for a star with a given central density, due 
to a further softening of the EoS. 



The efficiency with which QMC allows one to model hyperon contributions highlights 
the importance of baryon structure. The lack of a phase transition between hyperonic 

QMC and NJL quark-matter indicates the importance of both hyperon degrees of free- 
dom and dynamic chiral symmetry breaking, both of which must be treated with care. 
A transition is possible when one or both of these factors is neglected, but the physics 
must be removed by hand. 

In order to achieve a Gibbs phase transition between a baryon phase and a quark- 
matter phase in which dynamical chiral symmetry breaking gives rise to constituent- 
quark masses, a moderately stiff baryon EoS is required, otherwise no such transition is 
possible. Such is the case for hyperonic QMC and NJL quark-matter. 

While the QMC model at Hartree level successfully reproduces many properties of 
finite nuclei and stellar objects, the softness of the hyperonic/hybrid EoS (as evidenced 
by the maximum stellar mass) is unable to fully account for the most massive pulsars 
observed to date. We note that further improvements such as calculations at Hartree- 
Fock level (of particular interest are the Fock terms corresponding to the Tl meson) are 
currently in development. 



5. STELLAR SOLUTIONS 




(6) 



6. CONCLUSIONS 
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FIGURE 2. (Color online) Predicted total stellar mass and radius for hyperonic QMC with and without 
a phase transition to quark matter modelled with the MIT bag model Also shown are the data points from 
Refs. [8] (EXO 0748-676), [9] (EXO 1745-248), and [10] (4U 1608-52). 
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